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Abstract 

The three-dimensional half-shell t-matrix for a sharply cut-off Coulomb potential is analytically derived 
together with its asymptotic form without reference to partial wave expansion. The numerical solutions of 
the three-dimensional Lippmann-Schwinger equation for increasing cut-off radii provide half-shell t-matrices 
which are in quite a good agreement with the asymptotic values. 
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I. INTRODUCTION 



This is a continuation of a previous article [l| where the exact analytical three-dimensional wave 
function for a sharply cut-off Coulomb potential has been derived together with the corresponding 
scattering amplitude (the on-shell t-matrix). In the screening limit that scattering amplitude 
converges to a sum of two terms. One is the expected pure Coulomb scattering amplitude multiplied 
with the standard renormalisation factor (2pR)~ 2lT) ; the other one is new and includes angular 
dependent phases e ± 2l P Rsm 2 ; which oscillate without limit for infinite screening radius R. As has 
been conjectured in [2j that second term would disappear after integration over some angular 
intervals in the sense of a distribution. 



We are now interested in the screening limit of the corresponding half-shell t-matrix. The 

n 

pure Coulomb force result for that object is well known 3]. Its derivation goes back to work 
by [4]. Especially its discontinuous property at the on-shell point is of interest. In 



this 

property has been discussed based on a sharply cut-off Coulomb potential and using a partial 
wave decomposition. Like in our previous paper we felt that a direct three-dimensional approach 
avoids possibly open questions in that treatment related to the correct summation of the infinite 
number of partial wave components. (See QjEj], where the difficulties are spoken out leading in fact 
to incomplete results). We therefore study the half-shell t-matrix now based on the exact three- 
dimensional wave function for a sharply cut-off Coulomb potential and investigate its screening 
limit. The details of derivation are given in Section|TIl Numerical solutions of the three-dimensional 
Lippmann-Schwinger equation for different cut-off radii are compared with the asymptotic values 
in Section IIII1 We summarize and conclude in Section IIV1 

II. THE HALF-SHELL T-MATRIX 

For a a sharply cut-off repulsive Coulomb potential (for instance for two protons) 

V R (r) = @(R-r)- (I) 
r 

the exact three-dimensional wave function inside the potential range is given by 



with 



*W = Ae^Fi-ir,, 1, i(pr - p ■ r)) (2) 



A (2tt)I F(-in,l,2ipR) " (3) 



The normalisation corresponds to the choice of — ^ r e L ^' r as incoming wave. Further 77 = ^r— for 

(2tt)2 ' Zp 

two particle with mass m. Then the half-shell t-matrix is defined as 



<?'\Vr\*£ ) >=jA^ A J d\e~^- r y R (r)e^F(-i V ,l,i(pr-p-r)) . (4) 

We use the integral representation for the confluent hypergeometric function 

F(-i V , 1, i(pr -p-r)) = C(-i V , 1) / dt^—^r-e i ^ r -^ t (5) 

Jr t t 



with 

C(-iv, 1) = ^ (6) 

and r a closed path in the complex t-plane encircling t = and t = 1 in the positive sense. 
The r-integral is straightforward leading to 

with 



= y p 2 /; 2 - 2tp- A + A 2 , (8) 
A = p-p ' . (9) 



Thus 
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<#'|y R |4 +) >= -27r^^^C(-i77,l)y (10) 



and 



1 _ / 11 p i(pt+n)R _ 1 p i(pt-Q)R _ 1 



r 



Here we like to distinguish the two cases p' > p and p' < p and start with p' > p, where one 
can split (fTTj) as 



/• I _ f . I 1 p i(pt+n)R i(pt-Q)R r 1 + 1 1 

r = / r *<— )*'ln'wsr - + 2 h—r-tW^ (12) 



since the poles of ^tzqs do not lie between t = and t = 1. One has 

1 _ 1 

p2t 2 -fi2- 2 p-A(7i-4^) ' 



(13) 



It is easily seen that for p' > p the pole position to = as a function of p • p ' is always larger in 
magnitude than 1 if p > p'p ■ p' and smaller than zero if p < p'p ■ p' . Therefore we can choose from 
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FIG. 1: The original path of integration T used in Eq. (fT5|) . 



the very beginning the path T such that the pole at t = to nes outside that closed path. Since the 
only singularities of the integrand is the logarithmic cut between t = and t = 1 and the pole at 
t = we choose the T like in [1]. For the convenience of the reader this is depicted in Fig{TJ 

The second term in f)12|) is easily evaluated. If T is a circle with infinite radius the integral is 



zero. But changing the path T to that circle one picks up a residue due to (|13j) . If to : 
then 1 — to = [1 — to|e J7r and to = |to| ; if to < then 1 — to = |1 — to|e 2l7r and to = |to|e 
'^) ir, e-^. Consequently 

„/2 1 



v2 



> 1 



2p-A 

In both 



cases 



2 / (it(i_*)**i_ 

It t t p z t z 



1 



ft 2 



4vri b 2 -/ 
A 2 A 2 ' 



(14) 



Kow lu (he (irsl term in I 12), denoted as Y\. Again as in [jj we nplil I he integral between / f 
and t = 1 — e into two parts, choosing for instance t = 1/2 as intermediate border, and perform a 
partial integration for the integral J^ 2 to remove the pole |. In this way we get 

-l-e 



Y-l 



+ 

zero J one 



+ (1 



)( 



+ 



(15) 



and 



i _ t 11 P i(pt+n)R p i( P t-n)R 

dt{—^LL\l f 1 



t 7 t fr pt + ft pt-n 1 

_1 1_/ 1 p «(pi+n)_R p i(pt-U)R i 

_fi .y^Arf: _ ip 

i»7 v t ; n l pt + n pt-n J|e 

i /4 j i p i( P t+n)R p i(pt-n)R 
+ A. / dtt~ ir i—(l-t) iri —\- - 1 

+ i V J e dt { ' n [ P t + n pt-n 1 



(16) 



The integral f zero around t = is easily evaluated: 



dt( 



1 



f I -i i{pt+Q,)R 



zero t tQ pt + Q 

^cosARe- ir >-(l - e 
A z r] 



j(pt-n)R 
1 

pt-n J 



-27rm 



(17) 



It cancels exactly against the lower limit contribution at t = e in ([16 



The integral j vanishes as 0(e) and the upper limit (1 — e) in the last integral in (|15p can be 
replaced by 1. Thus as an intermediate result we have 



_ x 1 e i(f+n')-R e *(f-n') 



Y-\ = (I - e~ 2nv )(—) — \ 



2 2 

1/4 ^ .i e i(pt+n)R e i{pt~Q)R 



Q pt + Q pt-Q 
>i i _ / 11 i(pt+n)i? p i( P t-n)R 

+ (1 _ e -2^) / dt (:L_WAr£ _ £ i (i 8 ) 

v ' Ji y t ' tQpt + Q pt-Q 



2 



with 



f]' = y^_p.A + A 2 . (19) 

The lower integration limit e could be replaced by since only integrable logarithmic singularities 
remain. The differentiation in t leads to several parts. Only the pieces proportional to R will survive 
in the screening limit as will be argued below. They are given as 

P r\ 1 — t 1 p iR(pt+Q) JCt JR(pt-n) jq 

(1 - e_2 *" ) fi ^V'^WO-VrrO 1 ' (20) 

The contribution from the lower limit t = is evaluated by the method of steepest descent for 
R — ► oo . We use 

e iR(pt±U) _ e ±iRA e iRt(p T ^) ^ (21) 

P± > = "T' (22) 

and obtain 

- C 1 " e ~^)^ r ( 1 " ^V* A (P " + e" iRA (p + . (23) 

Now all contributions related to the arbitrary t = 1/2 border should cancel each other. This 
is indeed the case. By the same method of steepest descent the asymptotic contribution from the 
upper limit t = 1/2 in (|20p can be gained substituting t = 1/2 — r and using 

R *«CP*±n)< = ±ii?(n'±p/2) -j/?r(p=F^) 



^k=i/2 - jp , (25) 



as 



(i - e %y ^[wn-^+^-^rr^- 



1 _9™s int., e e 



^ 1 - e -'' v ^WTv-¥I^r h (26) 

This cancels exactly against the first part in (|18p . Further the last integral in (|18p contributes at 
the lower limit t = 1/2 for i? — > oo 

1 — t • 11 e*( p * +s ^ K — 1 e i{pt-0)R _ i 



;i -e~ 2 ^) J^dt^—^)^ 



tn l pt + vt pt-n 



It is easily seen that for |p • p'| ^ 1 the exponents p ± ^\t=i/2 7^ an d therefore that limit 
is O(i) and can be neglected. By analogous steps one finds that also the contribution from the 
upper limit t = 1 of that last integral in (|18p is 0(-^) in the screening limit. Therefore we obtain 
from (I18p and (|23p adding the remaining parts of the differentiation 

Y x -+ (1 - e" 2 ^)^r(l - ir 7 )i? i "[e i - RA (p - ^P + e" i/?A (p + 

ryA^ A A 

IT] Jq t u 

ir, Ida e *(pt+n)a e <(pt-n)Ji e i( P *+n)B d o eJ (pi-n)i? dn 

^ ( ~T^t ~ n dp pt + o pt - n J ~ (pt + o) 2 (i? + dp + p - n) 2 (p ~ dr J ' (28) 

Since there are no vanishing denominators nor pt =L = inside the range of integration the 
remaining integral is O(j^) and one ends up with the screening limit 

Y x -» (1 - e- 2m? )^r(l - iri)^\e iRA ip - + e" iRA (p + . (29) 

r)A z A A 

Then taking together with (|14p one finally arrives at 

y - < 1 -^ 2 *')^ r ( 1 -"')«"'[ e " ia (i'-^ : ) , " + ^" iA (p + ^)"] 

Thus like for the screening limit of the on-shell scattering amplitude given in l[ there result 
two terms, one, as expected, /^-independent and another still dependent on R. 

Now we turn to the case p' < p and start again from (jlip , which for a suitable path T can again 
be brought into the form (|12p . In this case the pole to from (|13p lies on the real axis between t = 
and t = 1. 

Since the path T encircles the cut the integral 
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FIG. 2: The modifications of T used in Eqs. f33jl -(|47 t . 



as is trivially seen by replacing the path T by a circle with infinite radius. 

Thus we are left with the first term in (|12j) which can be brought into the form 



Y ! dt( l - t Y r > 1 1 t ei(pt+n)R e i( ^-^ R (pt + fi) 
= ir ^ ^ pi + ft 2 p- A(t--^) 

^ v 2p- A y 



(32) 



We deform the path T into V such that the lower part of T is moved between and 1 into the 
upper half plane, as shown in Figj2j Thereby the path crosses the pole at to = ~^r, which leads 

2p- A 

to a residue: 



/ 

J pole 



1-t lv l 1 e^-Wjpt + n) Am p 2 -p' 2 iv 

y K 2p-A J 



(33) 



We used the fact that at the pole pt = ft. 

Next we further deform the path V such that it encircles t = coming from ioo and returning 
back to ioo in the positive sense and encircling t = 1 again coming from ioo and returning back to 
ioo in the positive sense. This new paths To and T\ are displayed in FigEl 

We get 



Y 



1 _ + 11 p i(pt+n)R i(pt-Q)R 
dt(- -) ir t--\- - 1 

ro v t ' tn [ pt + n pt-n J 



_ p 2 - p' jT? 

A 2 A 2 ; ' { } 

Because of the pole at t = we separate the integral J r into three parts (see Fig|2]) which read 

1 _ f 11 p i(pt+fl)R „i(pt-Q)R 

dt(- -) ir >-—\- - 1 

ro v t ' t n 1 pt + ft pt-n 1 

cft... + (l- e -*^) / dtt-^-\l-tr-{— — - — — ] . (35) 

Ji e ft pt + ft pt - ft 

The integral f zero is 

,l-t,„ll r e« f+!! ) fi ^( P t-n)R ,f e - 2m? j ]_ iAB e -iAR 

dt( — y V T^l-TT^ 3— FT = / . «e^#— r X -t- + 



JCT'tl 



2 



i 7 ift L pi + ft pt-n 1 J_3k (ee^yv ee^ A L A A 

A^ 77 



6-^(1 _ e-^^cosARe'^- . (36) 



3 1 J(pt+Q)R i{pt-Q)R 

(l-e" 2 ^) / dtt-il-^l-t)**- 



If 



Further 



ft pi + ft pt - ft 

1 1 i( P t+n)R i( P t-n)R 

= (i _ e -a» F)[ __ t -*, (1 _ -[__ -^-ir 

177 S Z pi + S Z pt — ft 

+ -/ _t)*J^[- — — ]] . (37) 

ir? 7 ie (ii v y ft L pi + ft pi -ft JJ v ; 

There is no contribution at the integration limit ioo and the lower limit contribution cancels 
against (f36j) . 

Thus we are left with the intermediate result 

A-rrj n 2 — n' 2 1 f io ° d 1 pi(pt+tyR i(pt-Q)R 



A-jri n — rf 1 P 00 d 1 f 



ft pi + ft pi - ft 

/■ i _ / 11 P i{pt+n)R i( P t-n)R 
The differentiation leads again to a piece explicitly proportional to 

p moo -i j. -i p i(pt+n)_R/ i t(pt-H)R/ rff?\ 



Y R = (1- e-*"»)- / *(—)"-[ TTTT^ 7 o ] ' (39) 

77 Jo * ft pi + ft pi — ft 

where we could put the lower limit of the integration to zero. The integral converges at the upper 
limit noting 

ft(ir) -» -irp+ ^— , (40) 

^ p • A . _ . 

dT^- p + V- (41) 
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Using 

e i( P t±n)R _^ e ±iRA e -RT( PT ££) ^ 
in the limit r — > we extract the leading behavior of Yr 

+ e -«^ e -a-(P+¥)(p + £^)] 

= (l_ e -*ni)i e -S? 1 [e*BAn(p_£l^) [ dTT -iV e -Rr(p-^) 
r] A 2 A y 

- (l-e- 2 ^)le-^i^r(l-iry)[ e « A (p-^r + e -^(p + ^r] . (43) 
Next we regard the integral from (|38|) 

/ dt(i — £)<*?£±r£ £ ] 

y ri v t ; tn 1 pt + n pt-n J 

/>l+ioo 1 _ / 11 p i(pt+U)R p i(pt—Q)R 

= <>--'•'>/ ^-'"iflW"^ 1 ' (44) 

To evaluate the leading contribution from the lower limit i = 1 we need 

e i(pt±n)R _^ e iR{p±p') e -Rr{pT^f) ^ 

and obtain 

f «... (1 _ e -^ )ie ^i_[^ [ dTe -^-^_^l f .(46) 
■M P 1 P + P Jo P-P Jo 

Since p ± = p(l ± p -p') ^ for ^ 0, 7r, which we exclude, the integrals in (|46p are 0(-^)- 
Finally the additional terms resulting from the differentiation in (|38p are given as 

-2™! r°^i-w r ,-^ 



(l-e-^)- / #(_-)«>-{(- , 
277 * s 2 1 — i i 2 ore 

e i( P t+n)ij e i( P t-n)R e i( P t+n)R e i( P t-n)R 

^ pt + n pt-n ^~ (pt + n) 2 ( - p + ^ + p - o) 2 ^ ~ It" ' ^ 47) 

It can be shown that along the imaginary i-axis the imaginary part of pt ± O is always positive. 
Therefore that integral, too, vanishes like O(^) in the screening limit. 
Thus we are finally left for p' < p with 

V _ 4 " p 2 - j/ 2 it> 

A 2 A 2 
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+ (l-e-^-e-^^R^^-ivW^ip-^^r + e- iRA (p+'^r] • (48) 
i] A 2 A A 

This is now to be compared with (I30p . repeated for the convenience of the reader and valid for 
p' > p 

Ani p' 2 -p 2 it? 
A 2 A 2 

+ (1 - e - 2 ^)-e§ r '-^ff ,? r(l - i V )[e iRA (p - + e- iRA {p + . (49) 

77 A^ A A 

We see that the ^-independent part jumps from p' < p to p' > p by a factor e~ 7r,? whereas the 
oscillating independent part jumps by a factor e 27rr? . 



Lastly we turn to the half-shell t-matrix given in (jlOp . The asymptotic value of A as derived in 

fl is 

A -» — ^-g- e^^pT^ + ifj) (50) 
(2tt)2 

and therefore the prefactor in (jlOp using ([6]) is asymptotically 

- 2ir-^jAC(-ir,, 1) - i^ef ^(2pi?)-^r(l + ir,) . (51) 
(2?r)2 (2vr)' 3 

This leads for p' < p to the half shell t-matrix element in the screening limit 



4 +) >-> ^e^( 2 ^)-^r(i + ifl )^i-j!L 



IT I 



- i£ri^^- P ^ + ^(^ + S^) r] , (52) 

where we used T(l + ir/)r(l - irj) = si ^ n ■ 

On the other hand the pure half shell t-matrix is well known [3j and given for p' < p by 

<p'Wc\*f +) >=^e^T(l + i V )±(^^)^ . (53) 
Therefore for p' < p we find the following result in the screening limit 

<p'\V R \^ ] >^e-^ ln2pR <p'\V c \^ {+) > 

2 1 [e m H\{p-^)r + e^ R \Up + ^)r\. (54) 



(2^) 2 A 2L v 2^ A 2 A 

The first term is the expected one as given in [5J. But there is, like for the on-shell t-matrix, 
an additional term, which only after integration over some angular region would disappear in the 
screening limit. 
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In case of p' > p the pure half shell t-matrix differs by a factor e nri and is 

< v'\Vc\*f +) >= pe"iT(l + iv)^(^^r ■ (55) 



Therefore in this case and using (|49p we find the following result in the screening limit 

- e 2 ^^^[e^ A (^-^)r + e-^ A (i(p+^)H . (56) 

The first term has the same structure as above, but the second one differs by the factor e 2wv 
from the one above. 

The first term has the same structure as above, but the second one differs by the factor e 27TV 
from the one above. 



III. NUMERICAL RESULTS 

It is interesting to compare the derived asymptotic forms (I54p and (I56p to the numerical solutions 
of the Lippmann-Schwinger equation for the sharply cut off Coulomb potential with different cut-off 
radii. As is well known [l|, Q] this equation can be written as a two-dimensional integral equation 

oo 1 

T(q',q,x') = ±-v(q',q,X f ,l) + J dq"q" 2 f dx"v{q\q\x\x")—^T{q\q,x"\ (57) 

o-i z ~ 

where 

„(„', ,y,„) = j4 v VW, q y X+ VT^VT^c 0Sv) (58) 



and m is the reduced mass of the system. 

For the sharply screened Coulomb potential of the range R considered in this paper 

rr/ i x e 2 1 - cos(QR) 

V(q ,Q,y) = ^2 Q2 > ( 59 ) 

where Q = \/q' 2 + q 2 — 2q' qy and the integral over ip in Eq. (I58p is carried out numerically. 

Solving the two-dimensional equation (|57p is a difficult numerical problem because V(q', q, x) 
shows a highly oscillatory behavior, especially for large R. We solved (|57p for positive energies 
where 

z = E c . m . +ie = ^ + ie (60) 
m 

11 



by generating the corresponding Neumann series and summing it up by Pade. In each iteration 
the Cauchy singularity was split into a principal-value integral (treated by subtraction) and a <5- 
function piece. All details about our numerical performance are given in PJ. By solving (|57p we 
obtain all matrix elements T(q',q,x; qo)', they can be chosen on-shell (as investigated in l|), half- 
shell or totally off-shell. Here we are interested in the half-shell elements, T{q' , qo, x; qo), and show 
examples in Figs. [3HH We choose just five (more or less arbitrary) values of x =-0.91, -0.50, 0.03, 
0.62 and 0.90, which corresponds to the following angles 9 between vectors q 1 and qo- & =155.5°, 
120.0°, 88.3°, 51.7° and 25.8°. Then for each fixed value of x we display the half-shell matrix 
elements T(q',qo,x; qo) as a function of q'. The analytical asymptotic forms of Eqs. (f54l) and (f56l) 
are given by the line and the numerical results are shown with symbols. We concentrate on the 
q' region in the vicinity of qo where the most interesting structures appear and skip the region of 
ligher q' values, where T{q' ,qo,x; qo) tends to zero showing more or less rapid oscillations. As in 
l| we restrict ourselves to the system of two protons scattering solely by the Coulomb force at 
E l p ab = 13 MeV. This gives q 0.396 fin -1 . Two cases of the cut-off radii R = 80 fm (Fig. (SJ and 
R = 200 fm (Fig. 0]) are considered. 

In the case of R = 80 fm the real part of T(q',qo,x; qo) is usually by one order of magnitude 
bigger than the imaginary part. The exception is x = 0.62, where the both parts are comparable. 
The analytical asymptotic form agrees rather well with the numerical result for the real part. The 
agreement is in fact very good in the region of q' < qo and a bit less satisfactory for q' > qo- For 
the imaginary part there are clear deviations between the analytical and numerical results, which 
become more pronounced for x > 0. In particular the analytical results show much more oscillatory 
behavior for q' > qo- Note also a sharp structure around q' = qo which develops for the imaginary 
part at x > 0. 

For R = 200 fm the real part of T{q' , qo, x; qo) is clearly dominant for all the considered values of 
x. The analytical asymptotic form shows much more oscillations than for R = 80 fm. It is clear that 
in order to trace these oscillations, many more points in the numerical solution would be required. 
Despite this fact, one can see at least fair agreement between the asymptotic analytical results and 
numerical solutions at the calculated points in the case of the real part. For the imaginary part, 
like in the case of R = 80 fm, the agreement is worse. This is presumably caused by limitations of 
our numerical treatment. 
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FIG. 3: (color online) The real (left) and imaginary (right) part of the half-shell t-matrix T(q', qo, x; qo) for 
the sharply cut-off Coulomb potential with the cut-off radius R — 80 fm. The solid (black) line represents 
the asymptotic analytical expression given in (|54|) and ([56]). The (red) dots show our numerical results. 
From top to bottom five different values of x are chosen: (a) x = —0.91, (b) x = —0.50, (c) x = 0.03, (d) 
x = 0.62, (e) x = 0.90. 
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FIG. 4: (color online) The same as in Fig. 3 but with the cut-off radius R — 200 fm. 



IV. SUMMARY AND CONCLUSIONS 



We investigated the screening limit of the exact analytical three-dimensional half-shell t-matrix 
for a sharply cut-off Coulomb potential. We used the exact three-dimensional wave function for 
a sharply cut-off Coulomb potential derived in Our direct three-dimensional approach avoids 
problems related to the summation of the infinite number of partial wave components. Numerical 
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solutions of the three-dimensional Lippmann-Schwinger equation for large cut-off radii agree fairly 
well with the asymptotic values. 
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